PROCEEDINGS 

OF SCIENCE 



43 




Bulk and finite-temperature transitions 

in SU(3) gauge theories with many light fermions 



(N 
O 

CN David Schaich*, Anqi Cheng, Anna Hasenfratz and Gregory Petropoulos 

Department of Physics, University of Colorado, Boulder, CO 80309 
Email: schaich@pizero.colorado.edu 



We investigate finite-temperature transitions in SU(3) lattice gauge theories with Nf = ^ and 
12 staggered fermions in the fundamental representation. For both of these systems, we have 
observed a strongly-coupled lattice phase in which the single-site shift symmetry of the staggered 
p) I action is spontaneously broken. Here we report new results for finite-temperature transitions 

on 24^ X 12 and 32^ x 16 lattice volumes, contrasting the 8- and 12-flavor systems. While the 
Nf = 12 finite-temperature transitions accumulate at the bulk transition bounding the strongly- 
^-H coupled lattice phase, the A^/ = 8 finite-temperature transitions are able to pass through the bulk 

^ transition, and behave as expected for a QCD-like system. We discuss our current results and the 

work in progress to complete our investigation of the finite-temperature phase diagram. 
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1. Introduction 

We recently reported the observation and characterization of a novel phase in SU(3) lattice 
gauge theories with A^/ = 8 and 12 nHYP-smeared staggered fermions [1]. From a variety of ob- 
servables (including the meson spectrum, static potential, low-lying eigenvalues of the massless 
staggered Dirac operator, renormalization group blocked plaquette and Polyakov loop, and newly- 
developed order parameters), we established that the single-site shift symmetry ("5^") of the stag- 
gered action is spontaneously broken (">S^") in this phase. In terms of continuum synmietries, the 

phase possesses both chiral symmetry and axial U(1)a symmetry in the chiral limit, even though 
the Polyakov loop and static potential clearly indicate confinement. However, we argued that the 

phase is likely to be a purely lattice phase with no continuum limit, on the grounds that 

1. its combination of confinement and chiral symmetry is forbidden by the continuum 't Hooft 
anomaly matching condition; 

2. it is bounded by first-order bulk (zero-temperature) phase transitions; and 

3. it appears in both 8- and 12-flavor systems, which we believe exhibit different infrared dy- 
namics (QCD-like spontaneous chiral symmetry breaking and IR conformality, respectively). 

Similar observations were also reported (with somewhat different interpretations) by other groups 
exploring the 12-flavor system with different staggered actions [2, 3, 4]. 

We have since extended our investigations of the 8- and 12-flavor systems through several 
complementary analyses, two of which (studies of the Dirac eigenmode scaling and Monte Carlo 
renormalization group) are discussed in other contributions to these proceedings [5, 6]. In this work 
we consider finite-temperature transitions in the 8- and 12-flavor systems, focusing on their behav- 
ior around the bulk transition that separates the >S^ phase from the weak-coupling phase connected 
to the continuum limit. Working primarily with lattice volumes xNt = 24^ x 12 and 32^ x 16, 
we observe an interesting contrast between the 8- and 12-flavor systems, shown in Fig. 1. For 
A// = 12, the finite-temperature transitions for A/^ = 12 and 16 are indistinguishable, and fall on top 
of the bulk transition bounding the ^phase. For Nf = 8, the finite-temperature transitions can pass 
through the bulk transition, moving to weaker couplings (larger jip) as Nt increases, in agreement 
with RG scaling. However, the 8-flavor system behaves differently at the lightest mass we consider, 
m 0.005, than at m > 0.01. We discuss this situation in Section 3. 

In Section 2 we present the observables that play the most prominent role in establishing the 
results shown in Fig. 1. First, the >S^ order parameters we introduced in Ref. [1] are the most robust 
means to resolve the bulk transitions around the JS^ phase, for T > as well as zero-temperature 
systems. Next, we review the validity of our procedure to improve the signal in the Polyakov loop 
by measuring it on RG-blocked lattices. Finally, the Dirac eigenvalue density p(A) is sensitive to 
the chiral properties of the system, providing additional information complementary to that from 
the other observables. In addition to distinguishing between chirally broken and chirally synmietric 
systems, we can also use p(A) to identify the phase. We conclude in Section 3 by discussing 
the implications of our current results, and our plans for completing this investigation. 

2. Observables for bulk and finite-temperature transitions 



The staggered action possesses an exact symmetry under the single-site shift transformation 
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Figure 1: Phase diagrams in the jip-m plane for A^/ = 8 (left) and 12 (right). The phase is shaded and 
the bulk transitions bounding it are shown in black (including preliminary 40^ x 20 results). The colored 
points indicate finite-temperature transitions determined primarily from the RG-blocked Polyakov loop and 
eigenvalue density. (The red diamonds at stronger couplings are determined from (xj/Xj/) .) 



where £,i^(n) = (— l)^v>Ai^v order parameters of this symmetry are 

AP^ = (ReTr D,,^ - ReTr □.+m,m>., even (2.2) 
AU^ ={a^{n)x{n)U^{n)x{n + ^) - a^{n + li)x{n + li)U^{n+ (2.3) 

In these expressions D^^^ indicates the plaquettes originating at lattice site n that include links in 
the jU direction; U and X are the gauge and fermion fields, respectively; OCi_i{n) = (— l)I^v<Ai"v. ^^yd 
the expectation value (• • •)^^ even taken only over sites whose /i component is even. 

In Ref. [1] we showed that AP^ and AUfj vanish in both the weak-coupling phase and the 
familiar chirally-broken lattice phase at strong coupling. Nonzero expectation values (in one or 
more directions /i) characterize the intermediate >S^ phase where the single- site shift symmetry is 
spontaneously broken. Even though the JS^ phase is bounded by bulk (zero-temperature) phase 
transitions, we observe the same transitions in the JS^ order parameters on the finite-temperature 
lattices we consider here. That is, 24^ x 12 calculations show the same transition in AP^ and A(7^ 
as do 12^ X 24 systems, and similarly for 32^ x 16 compared to 16^ x 32; these results are combined 
in Fig. 1, labelled "L = 12 /Nt = 12", etc. Our preliminary 40^ x20 data reveal only the transition, 
not yet the finite-temperature transition. 

In Fig. 2 we show the euclidean norms VAP 'AP and VAU-AU on 24^ x 12 and 32^ x 16 vol- 
umes for both the 8- and 12-flavor systems at fixed fermion mass m = 0.01. While these euclidean 
norms are not themselves order parameters, they are sensitive to the >S^ transition because they can 
only be large if AP^ and Ai7^ develop nonzero expectation values in one or more directions. The 

transition is volume-independent, with 24^ x 12 and 32^ x 16 results on top of each other. 
2.2 RG-blocked Polyakov loop 

Fig. 2 also includes RG-blocked Polyakov loop results, which show the A// = 8 finite-temperature 
deconfinement transitions moving as Nt increases, while those for A// = 12 do not. In pure-gauge 
systems, the Polyakov loop (TrL) is an order parameter, and it remains sensitive to deconfine- 
ment even in the presence of dynamical fermions. However, (TrL) becomes small and noisy as 
Nt increases, which motivated us to construct the improved observable (TrL/,) by measuring the 
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Figure 2: Signals of transitions in Pf with fixed m = 0.01 for A// = 8 (left) and 12 (right). The solid lines 
correspond to 32^ x 16 volumes, dotted lines to 24^ x 12. The green crosses and red squares are derived 
from the order parameters sensitive to the bulk transition bounding the phase. The blue triangles are 
RG-blocked Polyakov loop results, one indicator of the finite-temperature deconfinement transition. 

Polyakov loop on RG-blocked lattices. (TrL^) has the usual Z3 symmetry in the pure-gauge theory, 
and can also be thought of as an extended observable on the original, unblocked lattices. 

Our RG-blocking transformation consists of two HYP smearings with a = (0.6,0.2,0.2), as 
discussed in Ref. [7]. The once-blocked (TrLi) is therefore equivalent to the usual Polyakov loop 
measured with smeared links. As the number of blocking steps increases, the effective A^^ of 
the lattice is repeatedly halved, amplifying the signal in (TrL/,). To check that the RG-blocked 
Polyakov loop behaves appropriately, we measured it on Nf = 2+1 configurations made available 
by the HotQCD Collaboration [8]. The results in the left panel of Fig. 3 show that the RG-blocked 
Polyakov loop captures the same physics as the familiar observable, the only change being order- 
of-magnitude enhancements in the signal. At the QCD transition temperature Tc ^ 155 MeV, the 
maximally-blocked Polyakov loop on these 48^ x 12 lattices is (ReTrL2) = 0.142(1), while the 
unblocked value is (ReTrL) = 0.00188(13). 

The right panel of Fig. 3 shows maximally-blocked Polyakov loop data for the Nf = 12 system 
at pf = 2.7 and m = 0.01 on volumes 12^ x6 through 32^ x 16. Even though the unblocked Polyakov 
loop on these lattices is indistinguishable from zero for A/^ > 12, the RG-blocked data indicate that 
the system is deconfined through at least Nt = 16. It is important to note that even the maximally- 
blocked Polyakov loop can still indicate confinement: this is the case for the A// = 8 system with 
Nf = 16 at Pf = 4.7 and m = 0.01. Fig. 2 shows that this point is on the weak-coupling side of the 
phase, but not yet deconfined. 

2.3 Eigenvalue densities 

While investigating the 12-flavor phase in Ref. [1], we found the low-lying eigenvalues of 
the massless staggered Dirac operator to be a particularly striking observable. In the phase we 
observed a soft edge (a gap in the infinite-volume extrapolation of the eigenvalue density p(A)), 
indicating restoration of both chiral symmetry and axial U(1)a symmetry in the chiral limit. Fig. 4 
shows nearly indistinguishable p(A) on four zero-temperature volumes in the phase, which 
extrapolate to the soft edge Aq = 0.0175(5). While a soft edge is unusual, p(0) = also for chirally 
symmetric systems at high temperatures above the finite-temperature transition. Chiral symmetry 
breaking is characterized by p (0) 7^ 0. 
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Figure 3: Left: RG-blocked Polyakov loop results from A^y^ = 2+1 ensembles generated by the HotQCD 
Collaboration [8], on a log scale. The red circles are the unblocked observable, green triangles are once- 
blocked (i.e., smeared), and blue inverted triangles are twice-blocked. Right: Maximally -blocked Polyakov 
loop data for the Nf = 12 system at Pf = 2.7 and m = 0.01 on volumes 12^ x 6 through 32^ x 16. 
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Figure 4: The Nf = 12 eigenvalue density p(A) determined from the 200 lowest-lying eigenvalues on each 
of four volumes in the phase (Pf = 2.6, m = 0.005). 



The eigenvalue density p(A) therefore provides information about both the bulk and finite- 
temperature transitions. This is illustrated in Fig. 5, which presents results for the 8- and 12-flavor 
systems with m = 0.01 on 32^ x 16 volumes. In the phase (dotted lines), p (A) shows a steep slope 
for both systems, and a clear gap for A// = 12. As the coupling becomes weaker, A^/ = 8 moves 
from the >S^ phase to a chirally broken system with p(0) ^ 0; the transition to a chirally symmetric 
system only occurs at even weaker couplings around Pf ^ 4.85. For Nf = 12, in contrast, we move 
directly from the phase to a chirally symmetric system. In both cases, p(A) agrees with the 
behavior of the order parameters and RG-blocked Polyakov loop shown in Fig. 2. 

3. Discussion 

The selected observables reviewed in the previous section are those that we find most use- 
ful for identifying both the JS^ bulk transition (order parameters and eigenvalues) as well as the 
finite-temperature transitions (RG-blocked Polyakov loop and eigenvalues). The results of these 
investigations are shown in Fig. 1, which we now discuss in more detail. 

The most striking feature of Fig. 1 is the contrast between A/y^ = 8 (for which the finite- 
temperature transitions can pass through the bulk transition) and Nf = 12 (for which the finite- 
temperature transitions accumulate at the bulk transition). This behavior is clearly visible in both 
of Figs. 2 and 5 for the case of m = 0.01. Fig. 2 shows the >S^ order parameters falling to zero at the 
same Pf on 24^ x 12 and 32^ x 16 volumes. For Nf = 12 the RG-blocked Polyakov loop also jumps 
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Figure 5: Eigenvalue density p(A) for fixed m = 0.01 on 32^ x 16 volumes for Nf ■ 
(right). The dotted lines indicate [Jf in the phase. 
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to a large value at this Pf, while for A^/ = 8 (TrL^) only rises at weaker couplings. Similarly, in 
Fig. 5 we observe p (0) 7^ (chiral symmetry breaking) for 8-flavor systems on the weak-coupling 
side of the phase. The 12-flavor system on volumes up to 32^ x 16 is deconfined and chirally 
symmetric (p(0) = 0) for all couplings weaker than the ^ phase. 

A puzzling feature of Fig. 1 is the different behavior of the 8-flavor system at m = 0.005 
compared to larger masses. For m > 0.01, the change in Pf between the Nt = 12 and 16 transitions 
is Aj8/7 ^ 0.2, in rough agreement with the two-loop renormalization group prediction Aj8/7 ^ 0.25. 
At m = 0.005 the shift is much smaller (Aj8/7 ^ 0.05), and the transition in the order parameters 
moves the same amount instead of being volume-independent as at m > 0.01. 

In Fig. 6 we show the 8-flavor results at m = 0.005 for the order parameters, RG-blocked 
Polyakov loop and eigenvalue density, corresponding to that in Figs. 2 and 5 for m = 0.01. In 
contrast to both the 8- and 12-flavor results in Fig. 2, the transition in the JS^ order parameters 
for A/y^ = 8 at m = 0.005 moves slightly as Nt increases from 12 to 16. The rise in the RG-blocked 
Polyakov loop occurs at the same Pf where the order parameters vanish, unlike the A^/ = 8 panel 
of Fig. 2. Similarly, the m = 0.005 eigenvalue density in Fig. 6 shows no sign of chiral symmetry 
breaking for any J877; it appears to move straight from the phase to a chirally symmetric system 
withp(O) ^0. 

How should we interpret our A^/ = 8 results at m = 0.005? Fig. 1 shows that the 8-flavor finite- 
temperature transitions occur at steadily stronger couplings as the mass decreases. At Nt = 16, 
for example, the transition moves from Pf = 4.85 to 4.75 as the mass decreases from m = 0.015 
to 0.01. We suspect that at m = 0.005, the Nt = 16 finite-temperature transition is coincidentally 
located close enough to the JS^ phase that it becomes entangled with the transition separating the 
phase from the weak-coupling phase. This seems a more modest conjecture than the alternative 
that the nature of the 8-flavor system changes fundamentally for m < 0.005 compared to m > 0.01. 

We are currently generating 40^ x 20 gauge configurations to clarify this situation, which we 
will discuss in a future publication. If the Nt = 16 finite-temperature transition at m = 0.005 is 
indeed being affected by its proximity to the JS^ phase, then the corresponding Nt = 20 transition 
should be at weak enough coupling to scale in rough agreement with the renormalization group 
prediction for Aj8/7, similarly to A/^ = 16 at m > 0.01. We are also performing 40^ x20 calculations 
for the Nf = 12 system, to continue our program of directly comparing these two theories. 
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Figure 6: A^/ = 8 results for m = 0.005. Left: order parameter and RG-blocked Polyakov loop results as 
in Fig. 2. Right: eigenvalue density p(A) as in Fig. 5. 
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